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1. Introduction 



The Supersymmetrical Quantum Mechanics (SSQM) fj],[pf represents a concise algebraic 
form of spectral equivalence between different hamiltonian quantum systems realized by 
means of Darboux transformation 0,0. Certainly SSQM is used also for verification of 
some properties predicted by SUSY in the QFT |J. 

On the other hand the variety of extensions of SUSY-algebra seems to be broader in QM 
than in the QFT. Recently the polynomial SUSY-algebra has been proposed and elaborated 
in the one-dimensional SSQM 0. Its consequences for the scattering characteristics of 
spectrally-equivalent systems have been outlined in |7|]. 

The aim of our paper is to give the complete classification of second-order polynomial 
SSQM in one-dimensional and to generalize the polynomial SUSY onto two-dimensional QM. 
After short description of related superalgebra we select out the subclass of reducible SUSY 
algebras which is built by repetition of two single Darboux transformations or equivalently 
by gluing of two ordinary SUSY systems ||, (see also 0). We distinguish the irreducible 
subclass which cannot be factorized in this way. In the one-dimensional case both subclasses 
are rich of representatives as it is shown in Sect. 2. The two-dimensional version of the 
polynomial SUSY is built in Sect. 3 and it is found that the reducible subclass is restricted 
with potentials combining the harmonic oscillator and the centrifugal-type terms. 

The irreducible class appears to be more significant and is generated by the second deriva- 
tive supercharges. The set of nonlinear equations on potentials and coefficient functions of 
supercharges are derived. In Sect. 4 the corresponding binomial superalgebra is analyzed 
and nontrivial operator of dynamic central charge is discovered. When the senior deriva- 
tives of supercharge form the Laplacian the separation of variables is provided in particular 
coordinate systems. The explicit form of potentials is described and the operators of su- 
percharge and central charge are evaluated. In Sect. 5 the second derivative supercharges 
with metric gik(x) are introduced and their possible structure is analyzed. It is found that 
in the general case the polynomial superalgebra defines the dynamic symmetry operator of 
Hamiltonian system. Therefore the building of isospectral systems is tightly connected with 
another problem namely with a search for dynamic symmetry operator R. 



2. Polynomial SUSY in one dimension 



Let us remind the basic notations of the ordinary SUSY QM and describe its poly- 

nomial extension ||. The intertwining relations between two Hamiltonians with equivalent 
spectra are realized by means of well-known Darboux transformation [|3| : 

fcW = q+q- = -d 2 + V {1 '> (x) ; h {1 ~> (x) = E n ¥£ (x) ; 

h® = q-q + = -0 2 + V®(x); h^¥ 2 \x) = E n ¥ 2 \x)- 

h Wq+ = q+ h (*) ; q - h W = h (2) q -. (1) 

q + = —d + dx(x); q~ = (q + ) J< = d + dx(x); 
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where d = d/dx. 

The concise algebraic form of spectral equivalence is given by superalgebra for the part- 
ners hS 1 ' ,hy) and off-diagonal supercharges: 

{Q+, Q-} = H; {Q+f = {Q-f = 0; [H, Q±] = 0. (2) 

The superpotential x( x ) is defined by an arbitrary solution of the Schrodinger 

equation = E^(x) 

X{x) = -ln^(x). (3) 

If this solution is nodeless one has almost coinciding spectra of h^ 1 ' and hS 2 ' or, equivalently, 
double degeneracy of the energy spectrum of H. The polynomial superalgebra is created 
@ by the intertwining operators of higher order in derivatives. Let as consider the second- 
derivative superalgebra generated by the following operators 

q+ = d 2 -2f{x)d + b{x); q~ = (g+) f . (4) 

The conserving supercharges [H, Q ] = determine the Hamiltonian as follows : 

{Q + ,Q-} = (H + a) 2 + d, (5) 

where the potentials can be expressed in terms of a real coefficient function f(x) : 

= T 2f'(x) + f(x) 2 + Qr\-{ ££Vf - - a. (6) 



2f(x) \2f(x)J 4f(x 
The intertwining relations for H require that 

b(x) = -f'(x) + f(x) 2 ({PA* + j^. (7) 

2/0) \2f(x)J 4/(x) 2 

Depending on the sign of d one finds either the reducible algebra (d < 0) or the 
irreducible one (d > 0). In the reducible case there exists an intermediate Hamiltonian h 
which is a superpartner of both h^ 1 ' and with the ordinary superalgebra, respectively 
the second-order Darboux transformation can be factorized into a product of two ordinary 
Darboux transformation (see details in 0,0]): 

q + = qM = (-9 + W 1 )(-d + W 2 ); 

W 1 (x) = dxi(x); W 2 (x) = d X2 (x); (8) 
h W = qtqi + c/2; h = q^q+ + c/2 = q+q 2 - c/2; h {2) = q 2 q+ - c/2; 

which corresponds to a = 0, c 2 = —4c?. In the irreducible case (d > 0) c becomes imaginary 
and the hermitian intermediate Hamiltonian h does not exist. In fact the analytic continu- 
ation from c to ic leads to a complex potential for h. Respectively W\ and W2 are complex 
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functions and qf = =F<9 + Wi (i — 1,2), but qf ^ {q i )t. Nevertheless the second-order 
supercharges are hermit ian 

qfqt = (<fe (fif- (9) 
From the hermiticity of and Eqs.(^),(^) it follows that 

f'(x) r 
ReW 1 = -^tHt + f{x)\ ImWx 



2f(x) ' JV 7 ' 4/(1)' 

where /(#) is a real function as before and potentials are described by Eq.(Q). 

In general there is the large variety of pairs of potentials and which obey Eq.(§) 
for any sign of d. The higher-order polynomial SUSY algebra can be created by composition 
of ordinary (first-order in derivatives) SUSY transformations and of irreducible ones, which 
are of second-order in derivatives: 

{Q + ,Q-}= J] (H + c 3 ){{H + a t f + dj); d j > 0. (10) 

i+2j=n 

3. Two-dimensional Darboux transformations of 
second-order in derivatives 



The conventional two-dimensional SUSY algebra provides the equivalence of energy spec- 
tra between a pair of two scalar Hamiltonians hS 1 ' and 2x2 matrix Hamiltonian hik 
M = 1,2) & 

h w = _ d f + v m^) = -df + (d lX ) 2 + d?xi 

h ik = S ik df + S ik ((dix) 2 - dfx) + ZdidkX, 

where the summation on index I is implied, df = d\ + d\. 
The superalgebra is represented by intertwining relations: 

h (1) qf = qt^kh h ik q^ = q^h {1) ; qf = ^ + ^x(x) ; 
hikPk = Pjh (1) ] Pk h ki = h (1) pf, 
qfpf = 0; Pi = (pf) ] = e ik qt; e ik = -e ki . 

The energy spectra of h^ 1 ' and are different except for particular cases [IT|,[n . 
However it is of interest to extend this superalgebra so that the spectral equivalence could 
be established between some scalar hamiltonians. One can attempt to build the polynomial 
superalgebra similarly to the one-dimensional reducible scheme : 



h(1) ^ h fcW A T ( h(2) 



In result we arrive to the following intertwining relations between upper components 
and 

hWq + = q + h,W; q'h^ = fc<V, (11) 
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where 



h {2) = qt<li ~ 4d; q + = qfq^ = + ^x)(^< + $X)l 9 = 



(12) 



We have investigated the consistency conditions for superpotentials and have found the 
following solutions in this reducible case 



X 
X 



1 c 

2 ap2 ~ ~ ln P + ^^ 



1 2 C 

-ap 



V® = a 2 p 2 + 



Inp + + 6; 



r - e" + 



2 2 i 

a p + — 

P 2 



r 2 - e + 



2(a + c); 
2(a-c), 



where a, 6 are constants and is arbitrary function of angular variable <p. 

In two dimensions the irreducible class of second- derivative superalgebra with intertwin- 
ing components 



q 



-A + CA + B; q- = (q + y; A = df 



(13) 



can be constructed when imposing the conditions (ffl|) irrespectively of the existence of 
superpartners hS 1 ' and h^ 2 \ Thus we do not assume that q + = qfqi where [qtiQj] = 
[qf, q~j] = (the integrability condition). At this point the two-dimensional case is drastically 
different from the one-dimensional case. 

The intertwining between and leads to six equations on five real functions 
(C li2 ,B,VV>*>): 



)1U 



0; 



-d^C u - d v C„ + (V (2 - V^)S t 
AC V - 2d I/ (V {2) -B)- (V {1) - V {2) )C U = 0; 
-A(V {2) -B)- (V {1) - V {2) )B + C v d v V {2) = 0. 



(14) 



Certainly they cannot be satisfied for arbitrary potentials or V^ 2 \ Indeed the solu- 
tions can be found explicitly and read : 



C 



(C1 + 1C2) = az + 8/32 + 7; z = X\ + ix2] 

V {2) - B = \a I z I 2 +zp + z/3 + - I C \ 2 -77; 
4 4 

V&> - V w = d z C + dgC, 



(15) 



where a, r\ are real constants and j3, 7 are complex constants. The further evaluation of 
solutions for V^ 1 ' 2 ^ is based on the equation : 



(Cd z + Cd- Z )(B\C 



G I C 



where 



a 



G = a + {d z C){d- z C) ~t;(zC + zC) - 2(J3C + (3C). 



(16) 
(17) 



In terms of variables 



dz r dz . c dz r dz 
~c + 1 ~c ] lT2 = 1 ~C~ J ~C' 



one obtains 



C 



where Ffa) is an arbitrary function. Thus we find that the irreducible second-derivative 
SUSY algebra is realized in the class of potentials which is much broader comparing to the 
reducible case. 

4. The second-derivative superalgebra 

Let us describe superalgebra relations for the irreducible Darboux transformations (|TT|). 
Since supercharges Q ± (with components (jl3|)) commute with the Hamiltonian H one expects 
that the closing of superalgebra leads to the symmetry operator R (central charge) : 

{Q + ,Q-} = F(H,R); [H,R]=0. 

For the supercharges fll3|) one finds that 

F(H, R) = H 2 + R + 2rjH, (19) 

where R is a diagonal operator: 



Its components 



_ / R x 
R -[ R2 



R x = 2{a\z | 2 +A((3z + (3z))d z d- z - C 2 d 2 z - C 2 d\ - 

c(d z c)d z - c{d- z c)d- z + (Cd z + Od s )(B + v^) + 

B 2_ V (1)2_ 27]V (1). (2Q) 



R 2 = 2{a\z | 2 +A((3z + (3z))d z d s - C 2 d 2 z - C 2 d\ - 
C(d z C)d z - C{d- z C)d- z - (Cd z + Cd s )(B - V^) + 

B 2 _ y{2)2 _ 2T] y(2) (21) 

are the symmetry operators of Hamiltonians respectively. 

Let us derive the explicit solution of Eqs.(|l4|) and the related operators H, Q ± , R . They 
have different form depending on values of a and (3 in Eqs.(|l5D. 

i) If a = 0, (5 ^ 0, then 7 can be eliminated by translation z — > z — (7/8/?). In this case 
the variables 
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can be referred to the conventional parabolic coordinates [fT2 1 . From the Eq.([18|) one obtains 
that 

B = 2r l -\P\ 2 rf-2\f3\ 2 r 2 r 2 + F(r 2 ) 

and from Eqs.(^3j) one finds the following expressions for potentials 

y0)= -2n + l/?IV + FW (24) 



Tl + To 



vm= 2r 1±1 pp} + F M _ n (25) 

In terms of variables t±, t 2 the Laplacian is separable 

A = (d 2 + d 2 ) 

I P I 2 (t 2 + t 2 Y T1 T2 > 

with the same factor as in Eqs.(p^),(p5|). Hence the spectral problem for both Hamiltonians 
fity = Ety can be solved by i?-separation of variables (further on we omit the indexes 
of h,R,ty for brevity). Namely one can decompose the corresponding eigenf unctions for a 
given E into a sum 

^ = J2 ^n01n(ri)02n(r 2 ), (26) 

n 

where v n are constants and <f>i n (Ti), 2n (T 2 ) are solutions of one- dimensional equations : 

-K(n)+ I P I 2 HE + 17)7? =F 2r 1+ | /? | 2 t?] fc„fa) = ^i„(n); (27) 

I I' h(J5 + V)r 2 2 + F(r 2 )]0 2n = -^<Mt 2 ). (28) 

where the upper (lower) sign in Eq . ( 27 ) corresponds to h^ l >{hS 2 ' and \ n is a separation con- 
stant which serves as a spectral parameter for the symmetry operator R : -R0i n (Ti)0 2ri (r 2 ) = 
(A n — 77 2 )0in(Ti)02n(' r 2)- Depending on properties of the function F(r 2 ) one can obtain ei- 
ther the finite-dimensional subspace {Aj} of degeneracy for the energy spectrum or even the 
infinite-dimensional one. However the spectral values of h and R form the complete set of 
quantum numbers (the integral of motion) which characterize uniquely the quantum state. 

ii) Let us describe the second case (3 = 0, a > when the appropriate coordinates are 
elliptic ones: 



n = -±=ln(z + Jz 2 + 1 )(z + Jz 2 + 1 ) ; (29) 



a \ v ay \ V a 



i z + Jz 2 + * 

r 2 = ~L/n (30) 

Va z + ,/z 2 + ^ 



In a full analogy to case i) one can find 



5 = ;«T7y (29 " ;i " \ sl ~ hh + F(T2)) '' (31) 
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,/(1, = 5aTTA) ( - 29 " /l + 5 /f + F(T2)) -" ; (32) 

V' 2> = 2 (h + h) {2dr ' h + I 1 * + F(T2)) " 1- (33) 



where 



1 / I 7 I 2 \ 

fx = - f a ■ exp{^/ar 1 ) + — — exp{-^far x )\ ; (34) 

/2 = ^ (j ■ exp(iyfaT 2 ) + 7 • exp(-iy/ar 2 )^ . (35) 

This coordinates allow to separate variables as well that leads to the solution for the wave 
function in the form (EBl) where now 



- K(ri) + \ (rdrji + \fl -{E + v)f^Mn) = ^0m(ri); (36) 

-02n(r 2 ) + \ (p^- -(E + ij)/ 2 Wfa) = -^0 2 „(r 2 ) (37) 

and A n is again eigenvalue of the symmetry operator Eqs.([20l),(|21~l). 

The analysis of the case (3 = 0, a < is similar and its result can be formulated as 
follows. One has to replace real (ti,t 2 ) by imaginary (— zr 2 , — iT\) in Eqs.(p9|), (j30|) and 
respectively (/ifa), / 2 (t 2 )) by -/ 2 (r 2 ), -/ifa) in Eqs.(^), Then the relations for 

potentials B,Vi,V 2 are given again by Eqs.(|3~l~D - (|33|). 

We remark that the supercharge is factorizable into a composition of two ordinary SS- 
operators (qfq[) for the special choice a > 0, /3 = 7 = that corresponds to the separation 
of variables in the polar coordinates (see Sect. 3). We conclude that from the supersymmetry 
or equivalently from the intertwining of two Hamiltonians one inevitably recovers the hidden 
dynamical symmetry realized by R and furthermore the .R-separation of variables. 



5. General case and discussion 



Let us describe the natural extension of the second derivative Darboux transformation 
which is generated by the operator with metric 

q + = gik{x)did k + Cidi + B. (38) 

The intertwining relations (|T1|) determined completely the form of metric gik{x) which sat- 
isfies the following equation: 

digik{x) + digi k (x) + d k g u {x) = 0. 

Its solution can be easily found 

011 = ay 2 + aty + h; 
g 22 = ax 2 + a 2 x + b 2 ; 

022 = -~{2axy + a x x + a 2 y) + 63. 
8 



Thus one can see that in senior derivatives the possible supercharges belong to the E(2)- 
universal enveloping algebra |TJ| in which, at the level of second order, we distinguish three 
different possibilities : 



aP 2 + 7 A + CA + B; 

a J 2 + (3P 2 + 7 A + Qdi + B, 

where J and P are rotation and translation generators, respectively,and a^O. 

The coefficients of supercharge Eq.(p8|) and the potentials V^'^ 2 ' obey the modified 
equations: 

d t C k + d k Ct + Ag lk - (V w - V {2) )g lk = 0; 
ACj + 2d t B + 2g lk d k V^ - (yd) - V ®)Ci = 0; 
AB + g lk d k d t V (2) + CAiV {2) - (y« - V^)B = 0. 

The generalized superalgebra evidently yields to the symmetry operator for the Hamiltonian 
H: 

{Q + ,Q-} = R; [R,H]=0. 

However this operator creates the dynamical symmetry of higher order which cannot be 
in general represented by a polynomial of the Hamiltonian and of second order symmetry 
operator (similar to Eq. (|T9"D ) . But in any case the closing of SUSY algebra leads to the 
.R-separation of variables in the spectral problem for the Hamiltonian H (to the integrability 
of the corresponding dynamical system). 
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